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The low energy excitations of graphene can be described by a massless Dirac equation in two 
spacial dimensions. Curved graphene is proposed to be described by coupling the Dirac equation to 
the corresponding curved space. This covariant formalism gives rise to an effective hamiltonian with 
various extra terms. Some of them can be put in direct correspondence with more standard tight 
binding or elasticity models while others are more difficult to grasp in standard condensed matter 
approaches. We discuss this issue, propose models for singular and regular curvature and describe 
^^ ■ the physical consequences of the various proposals. 
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I. INTRODUCTION 



Since its experimental realization graphene has been a focus of intense research activity both theoretical and 
experimentally as can be seen in the recent reviews [l|, 0, S Ij]- The origin of this interest lies partially on the 
^ ', experimental capability to tailoring the samples into special geometries leading to graphene based electronics but it 
5-H is its unusual electronic properties and the breakdown of the standard Fermi liquid description what has lead the 

c/3 ^ main activity in the field. Despite the intense research three aspects of the electronic properties remain to be fully 
understood: that of the minimal conductivity p, [g, Q j the observed charge inhomogeneities at low densities [8|] , and 
2 ' the exceptionally high mobility at room temperature with an associated mean free path of the order of the size of the 
^ . samples j9i, ilfl] . A key point to understand all the three problems is a better knowledge of the nature and effect of 

' ' disorder in graphene. 

Zh . One of the most intriguing properties of the suspended graphene samples is the observation of mesoscopic cor- 

C^ rugations in both suspended 11, 12] and deposited on a substrate Il3l. 1 14l | samples whose possible influence on the 
electronic properties only now starts to be realized [ll|, [I^ Ua, [l3, \lS, \1m ■ Although the observed ripples were invoked 
from the very beginning to explain the absence of weak localization in the samples [20, l21| , there have been so far 
few attempts to model the corrugations. The two main approaches are based either on the presence of disclinations 
and other topological defects in the graphene lattice [22, [23| or on the theory of elasticity [llj, [l^, [la] ■ The concrete 
^\^ ■ realization of both models in the continuum gives rise to the appearance of gauge fields coupled to the electronic 
f^ , degrees of freedom [24| . In this work we will describe two different approaches to model curvature in graphene and 
0^ ■ study its influence on the physical properties of the material. 

cn ; 

O ■ II. A SUMMARY OF GRAPHENE FEATURES 

oo : 

^7^ ' Under a theoretical point of view the synthesis of graphene has opened a new world where ideas from different 

I> . branches of physics can be confronted and tested in the laboratory. On the electronic point of view it can be shown 

I^ ' that the low energy excitations of the neutral system obey a massless Dirac equation in two dimensions. This special 

^ ^N , behavior originates on the geometry and topology of the honeycomb lattice and has profound implications to the 

d ' transport and optical properties. Although the Fermi velocity is approximately a hundredth of the speed of light, the 

masslessness of the quasiparticles brings the physics to the domain of relativistic quantum mechanics where phenomena 

like the Klein paradox or the Zitterbewegung [2] can be explored. None of these questions arise within the quantum 

field theory approach but its full applicability to the condensed matter system is questionable. Monolayer graphite 

- graphene- consists of a planar honeycomb lattice of carbon atoms shown at the left hand side of Fig. [TJ In the 

graphene structure the in-plane a bonds are fornred fronr 2s, 2px and 2py orbitals hybridized in a sp^ configuration, 

while the 2pz orbital, perpendicular to the layer, builds up covalent bonds, similar to the ones in the benzene nrolecule. 

The a bonds give rigidity to the structure and the n bonds give rise to the valence and conduction bands. A usual 

tight binding analysis of the system [2a, [2al leads to a special band structure shown at the right hand side of Fig. [T] 

The Fermi surface of the neutral system consists of six Fermi points (only two are independent) . A continuum model 

for the low energy excitations around the Fermi points (i=l,2) can be defined governed by the Hamiltonian 



no^ = ihvF /"cf2r§,(r)[(-l)V,9, + aydy]^^{r) 



(1) 





FIG. 1: Left: The honeycomb lattice is made of two interpenetrating triangular lattices. Right: the dispersion relation. 



where a^^y are the PauU matrices, up = (3to)/2, and a = \AA is the distance between nearest carbon atoms. The 
components of the two-dinrensional wavefunction: 

correspond to the amplitude of the wave function in each of the two sublattices (A and B) which build up the 
honeycomb structure and the subindex i refers to the two Fermi points. It can be shown that the wave function ([2]) 
transforms as a bispinor in k-space and it acquires a phase of tt under a 27r rotation. This Berry phase together with 
the helicity conservation has important physical consequences as the absence of backward scattering in the electronic 
transport what implies the absence of weak localization. 

The spinor structures around each of the two Fermi points remain degenerate and independent in the absence 
of short range interactions of disorder. As we will see, topological disorder mixes the two Fermi points and in the 
modelling of them it will be useful to combine the two bispinors into a four dimensional object described by the 
effective Hamiltonian 

Hd = -ivFTT'{'^®(^idx + T^ ®(^2dy), (3) 

where a and t matrices are Pauli matrices acting on the sublattice and valley degree of freedom respectively. 

III. MODELLING CURVATURE IN GRAPHENE 

When trying to take into account the curvature of the graphene samples and its physical implications there are 
two aspects to consider: The first and of the most interesting questions is the physical origin of the ripples and their 
dynamics. The present experimental situation seems to indicate that there is little dependence of the corrugations 
with temperature -although systematic studies have not been performed yet -. On the other hand, bilayer structures 
also present corrugations but less pronounced [l5|. These aspects concern the sigma bonds of graphene and their 
elastic properties and involve energies of the order of tens of eV. Although the elastic and mechanical properties of 
carbon nanotubes have been extensively explored [27| , very little can yet be found on the elastic properties of plain 
graphene |28t . 

A second aspect is that of the influence of the ripples on the electronic properties of the samples. A possible 
approach to this problem is to assume that the samples are corrugated for whatever reason and devise a model to 
study the implications of the corrugations on the electronic properties. This is a sensible procedure considering that 
the electronic properties of graphene are attached to the tt bonds and the associated processes involve energies several 
orders of magnitude smaller than those related to the elasticity of the a bonds. 

Since the low energy excitations of graphene are well described by the massless Dirac equation, a natural way to 
incorporate the effect of the observed corrugations at low energies is couple the Dirac equation to the given curved 
background. The main assumption of this approach is that the elastic properties of the samples -determined by the 
sigma bonds - are decoupled from the (pi) electron dynamics. The ripples can then be modelled by a fixed metric 
space defined phenomenologically from the observed corrugations and the electronic properties of the system will 
be found from the computation the Green's function in the curved space following the standard formalism set in 
gravitational physics 29, 30]. 

Although the physical origin of the ripples is unknown, we could distinguish two possibilities: those coming from the 
curvature in the substrate [14] that can be modelled by a smooth curved space and those observed in the free standing 



samples. The last class needs to be modelled by including topological defects [3l|: disclinations and disclination 
dipoles (dislocations). In what follows we will summarize the general formalism to include a curved background and 
to extract the physical implications. 

The dynamics of a massless Dirac spinor in a curved spacetime is governed by the modified Dirac equation: 

z7''(r)V,,V = (4) 

The curved space 7 matrices depend on the point of the space and can be computed from the anticommutation 
relations 

{7''(r),7'^(r)} = 25^''(r). 

The covariant derivative operator is defined as 

V;, = 5^ - il^ 

where f2^ is the spin connection of the spinor field that can be calculated using the tetrad formalism [23] ■ 

Once the metric of the curved space is known there is a standard procedure to get the geometric factors that enter 
into the Dirac equation. In the modelling of the graphene ripples, the metric can be treated as a smooth perturbation 
of the flat surface and physical results are obtained by a kind of perturbation theory. Very often, the final result can 
be cast in the form of the fiat Dirac problem in the presence of an effective potential induced by the curvature. 

The electronic properties of the curved sample can be extracted from the two point Green's function. The equation 
for the exact propagator in the curved space-time is: 

ire^^'id. + n,) G{x,x')=S{x~x'){~g)-i, (5) 

where i7^ is the spin connection and \/—g is the determinant of the metric. From the propagator one can compute 
the local density of states: 

piE, r) = ^-ImTr[G{E, r, r)70], (6) 

TT 

and other one-particle properties of the system as the electron lifetime. 

The described formalism can in principle be worked out for any curved background. The practical difficulties to 
implement this scheme are related with the choice of a realistic metric parametrizing the graphene sheet and to the 
ability to perform a "weak field" expansion of the metric around the fiat case. In the next sections we will see the 
mechanism at work in two quite general examples. 

IV. SMOOTH RIPPLES FROM THE SUBSTRATE 

The case of smooth curvature was studied in reference [l7l] where a general metric was considered that was non- 
singular over the surface and asymptotically fiat. The embedding of a two-dimensional surface with polar symmetry 
-for simplicity- in three-dimensional space is described in cylindrical coordinates by a function z(r) giving the height 
with respect to the flat surface z=0, and parametrized by the polar coordinates of its projection onto the z=0 plane. 
The metric for this surface is readily obtained by computing 

dz'^i^fdr^^afir)dr\ (7) 

and substituting for the line element: 

ds^ = dr^ + r^dO^ + dz^ = (1 + af{r)) dr^ + r^dO^. (8) 

where /(r) is a smooth function with the appropriate asymptotic behavior and the parameter a controls the deviations 
from flat space. The advantage of this approach is that it provides already the perturbative parameter needed to 
compute the electron propagator. 

A concrete example is provided by the gaussian shape shown at the left hand side of flg. [21 It models a smooth 
protuberance fitting without singularities in the average fiat graphene sheet. 

/(r)=4(r/6)2exp(-2rV6'), (9) 





FIG. 2: Left: A smooth curved bump in the graphene sheet. Right: Effect of the curved bump shown in the left on the local 
density of states. The color code is indicated in the figure. Darker (lighter) areas represent negative (positive) corrections to 
the density of states of the flat graphene sheet. 



This example was worked out in full detail in ref . [17| . A comparison of the Dirac Hamiltonian in the plane (flat) in 
polar coordinates 



Hflat = hvF 



with the corresponding curved hamiltonian 





dr-l^ + i- 

I r 2r 



^curved — llVp 



dr + i^ + i- 

' r 2r 



{l + afir))-^/^dr^t^ + Ag 



{l + af{r))-^/^dr + i^ + Ag 



(10) 



(11) 



Ao = 



2r 



l-(l + a/)-l/2 



2r 



(12) 



allows to extract the two main features of the model. First we can see that the curved bump produces an effective 
Fermi velocity Vr in the radial direction given by 



«,(r,0)=z;f.(l + a/(r))-i/2. 

The effective Fermi velocity will always be smaller in magnitude than the flat one. 
described in polar coordinates by an arbitrary function z = z(r) it will be 

Wo 



v/l + ^'(r)2 



(13) 
For a general curved surface 



(14) 



In a more general case we will have the two components of the velocity changed but always to a smaller value 
irrespective of the sign of the curvature at the given point. 

The second feature that arises is an effective magnetic field perpendicular to the graphene sheet given by 



B, = --drirAg 
r 



1 



«/' 



4r (1 + a/)3/2 ■ 



(15) 



The electronic properties of the curved sample were computed from the electron propagator to first order in the 
small parameter a that measures the deviation from flat space. In our example, a — {A/b)'^ is the (squared) height 
to length ratio of the gaussian, so for typical ripples in graphene a « 0.01, since this ratio is of the order of 0.1 [1]| . 
They are shown in the right hand side of Fig. [21 As a general feature the curvature induces inhomogeneities in the 
electronic density of states that could be related to the observations reported in f8l|. It is worth noticing that the 
contribution of the effective gauge field coining from the spin connection to first order in perturbation theory vanishes 
and all the corrections come from the determinant of the metric and the curved gamma matrices. This is important 
when trying to compare the present formalism with the tight binding modelling of curvature where only the gauge 
field arises from the modulation of the hopping parameter. 



V. TOPOLOGICAL DEFECTS 

In the absence of a substrate or strain fields the only way to have intrinsic curvature in the samples is by the 
presence of topological defects [3i| . In the case of the hexagonal lattice of graphene disclinations form by replacing a 
hexagon of the lattice by a n-sided polygon with n 7^ 6. The most common disclination in graphene is the pentagon 
that plays an important role in the formation of fullerenes [S^- Dislocations are made of pentagon-heptagon pairs 
and they have been widely studied in connection with the properties of carbon nanotubes [27[ and, more recently, in 
graphene [33, [SJ, S, [SO] . Observation of topological defects in graphene have been reported in [33, [Sa, [33] . 

Following the discovery of the fullerenes and nanotubes topological defects have been modelled in the hexagonal 
lattice in different ways. A very interesting approach related with the theory of elasticity is the gauge theory of 
the defects described in refs. [40, El, 42., ji3\ or the metric formulation of the theory of defects in solids set in [44| . 
Most of these approaches model the topological defects focussing on their "holonomy" [43, |4a, [43, [43, |43]and very 
little effort has been devoted to the specific issue of the curvature in these cases. The curvature of the graphene 
surface in the presence of a disclination in the continuum limit has a delta-function singularity at the position of the 
defective ring. In modelling spherical [4^ |4^ and quasi-spherical [53, ImI fullerenes curvature effects are included but 
the curvature is averaged along the spherical surface. The more difficult task of modelling several topological defects 
located at arbitrary positions was undertaken in refs. [2^, |23|. An equal number of pentagon and heptagon defects 
was considered to keep the samples flat in average. The proposed metric was taken from the cosmic string scenario: 

where 

N 



A(r) = Y^ 4(1, login 



and 

n = [{x - a,f + {y - b,f]'/\ 

This metric describes the space-time around N parallel cosmic strings, located at the points {ai,bi). The parameters 
/ij are related to the angle defect or surplus ct by the relation Ci — 1 — 4/Xi in such manner that if q < 1(> 1) then 
fii > 0(< 0). Within the formahsm described earlier it can be seen that the Green's function is modified by the 
potential 

F(w, r) = 2iA-f°do + lA-f^dj + ^-i\djK). (17) 

The computation of the corrections to the local density of states for different positions of the defects showed that 
pentagonal (heptagonal) rings enhance (depress) the electron density, a result that was obtained previously [52] with 
numerical simulations. Numerical ab initio calculations show sharp resonant peaks in the LDOS at the tip apex of 
nanocones [S^, [53] that have been proposed for electronic applications in field emission devices. Similar results have 
also been obtained analytically in [54| . A recent calculation of the minimal conductivity in graphene with a random 
distribution of heptagon and pentagon rings bases on this model has been performed in [55| . 

VI. CONCLUSIONS AND FUTURE 

From the present and similar analyses we can conclude that the morphology of the graphene samples is correlated 
with the electronic properties. In particular the presence of ripples, irrespective of their origin, induces corrections to 
the density of states and affects the transport properties. The spatial variation of the Fermi velocity is a distinctive 
prediction of the given geometric formalism. We note that the effective Fermi velocity is the fitting parameter used in 
most experiments [8| whose interpretation might change if the possibility of a space-dependent velocity is considered. 

Gauge fields are abundant in graphene and arise in many different contexts. This is partially a proof of the 
robustness of the Dirac formulation where only minimal coupling in the form of electronic current-vector field are 
marginal interactions to be considered at low energies [5 61 ] . They arise in various contexts and it would be interesting 
to have a complete classification. In the present context it is interesting to note the work of [53] where it was shown 
that, around each Fermi point, static strains can mimic the effects of external electric and magnetic fields. 



Future work should address a better understanding of the mechanisms leading to ripple formation and elastic 
properties of graphene. There is also little work done on transport and localization properties of graphene in the 
presence of topological defects, a very important issue. 
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